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Abstract 

We calculate the baryon octet masses to fourth order in chiral perturbation theory employing dimensional 
and cut-off regularization. We analyze the pion and kaon mass dependences of the baryon masses based on 
the MILC data. We show that chiral perturbation theory gives stable chiral extrapolation functions for pion 
(kaon) masses below 550 (600) MeV. The pion-nucleon sigma term in SU(3) is also investigated, we find 
cr^]v(0) = 50.7... 53.7 MeV. 
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1 Introduction and summary 



The masses of the ground state baryon octet are of fundamental importance in the investigation of three-quark 
states in QCD. With the advent of improved techniques in lattice QCD and systematic studies within the 
framework of chiral perturbation theory, one can hope to gain an understanding of these quantities from first 
principles. Present day lattice calculations are done at unphysical quark masses above the physical values, 
therefore chiral extrapolations arc needed to connect lattice results with the physical world, provided that the 
masses are not too high (for an early approach to this problem, see e.g. pP). With the recent data from the 
MILC collaboration |21 El it appears to be possible to apply chiral extrapolation functions derived from chiral 
perturbation theory (CHPT) (for a review, see j^). In this paper, we analyze the baryon masses as a function 
of the pion and the kaon masses in CHPT based on the MILC data. Other pertinent lattice papers are e.g. 
El 13 IHl El- It should be said from the beginning that we ignore the efi'ects of the a) the finite volume, b) 
the finite lattice spacing and c) the staggered approximation in this study"^^ since our aim is more modest - we 
want to find out whether these chiral extrapolations can be used for the presently available lattice data. Once 
this test is performed, one should then apply the full formalism including the abovementioned effects. Note 
that the quark mass expansion of QCD is turned into an expansion in Goldstone boson (GB) masses in CHPT 
- we thus use both terms synonymously. 

The baryon masses have been analyzed in various versions of baryon CHPT to third and fourth order, for an 
incomplete list of references see ^1 ^1 ^1 El El El El ■ Our investigation extends the work in the two-flavor 
sector presented in j^U] and we heavily borrow from the earlier SU(3) calculations of El El- The pertinent 

results of this investigation can be summarized as follows: 

1) We have calculated the baryon masses to third and fourth order in the chiral expansion making use of 
cut-off regularization as proposed in |20| . As in that paper, we have also considered an improvement term 
at third order to cancel the leading cut-off dependence in the baryon masses, see Sect. 13.21 

2) The improvement term consists of three independent terms, whose cut-off independent coefficients have 
been determined by considering the nucleon mass (to allow for a direct comparison with the SU(2) calcu- 
lation of Ref. "W). We have demanded that for the physical pion and the physical kaon mass the nucleon 
mass passes through its physical value. This fixes two parameter combinations. The third parameter is 
determined from a best fit to the trend of the earlier MILC data [2] for m.M{M.,^), cf. Fig. El under the 
condition that the deviation from the earlier determination of the corresponding low-energy constants in 
El is of natural size. We find indeed a visible improvement in the description of the lattice data and also 
much better stability under variations of the cut-off. 

3) The full fourth-order calculation utilizing the low-energy constants as determined from the improvement 
term leads to an accurate description of the MILC data for pion masses below 550 MeV, see again Fig. El 
Note that the two lowest mass points of the more recent MILC data can not be well described. Also, 
the use of the low-energy constants (LECs) from El l^ad to a less satisfactory description. We have also 
discussed the theoretical uncertainty of this procedure, cf. Fig. El 

4) From the pion mass dependence of the nucleon mass, we can deduce the pion-nucleon sigma term. For 
the best sets of low-energy constants, we find crTrAf(O) = 50.7. . .53.7MeV. 

5) The kaon mass dependence of the nucleon mass is less well determined. Still, the extrapolation functions 
can be applied to kaon masses below ~ 600 MeV. We deduce that the baryon octet mass in the chiral limit 
lies in the interval 710 MeV < toq ^ 1070 MeV. This is consistent with earlier estimates, see e.g. {15] . 

6) We have also considered the pion and kaon mass dependences of the A, the E and the S and compared 
to the existing MILC data, cf. Figs. I5I6I7I Note that we have not fitted to these masses. Our chiral 
extrapolations for the E and in particular for the S as a function of the pion mass are flatter than the 
MILC data. This is partly due to our strategy of flxing all parameters on the nucleon mass. We remark, 
however, that one should expect a decreased pion mass dependence as the number of strange valence 
quarks increases. 

Some of these effects are studied in 1101 1111 1121 (and references therein) . 
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The material in this paper is organized as foUows. Section |21 contains the effective Lagrangian and a short 
discussion about the various regularization methods employed in our calculations. In section |31 the ground 
state baryon octet masses are given at third, third improved and fourth order in the chiral expansion. In 
particular, we concentrate on the differences to the SU(2) case [201 (we refer to that paper for many details). 
Our results for the various baryon masses as functions of the pion and the kaon masses and the stability of 
these results under cut-off variations are given and discussed in section 01 Many technicalities are relegated to 
the appendices. 



2 Formalism I: Generalities 

In this section, we display the effective Lagrangian underlying our calculations and discuss briefly the cut-off 
regularization utilized and its relation to the more standard dimensional regularization (DR) . We borrow heavily 
from the work presented in Refs. [T51l2U| and refer the reader for more details to these papers. 



2.1 Effective Lagrangian 

Our calculations are based on an effective chiral meson-baryon Lagrangian in the presence of external sources 
(like e.g. photons) supplemented by a power counting in terms of quark (meson) masses and small external 
momenta. Its generic form consists of a string of terms with increasing chiral dimension, 

r - -X- , /-(a) , /^(4) , /^(2) , /.(4) , . 

Here, B collects the baryon octet and (p stands for the Goldstone boson octet. The superscript denotes the 
power in the genuine small parameter q (denoting Goldstone boson masses and/or external momenta). The 
explicit representations of 4> and B are: 

^TT- -^° + ^7y x/2i^" , B{x)^\ E- -isO + ^A n I .(2.2) 



A complete one-loop (fourth order) calculation must include all tree level graphs with insertions from all terms 
given in Eq. (|2.1|l and loop graphs with at most one insertion from C'^^ or £^ . Throughout, we employ the 
heavy baryon approach, which allows for a consistent power counting since the large mass scale (the baryon 
mass ms) is transformed from the propagator into a string of 1/ms suppressed interaction vertices. The lowest 
order (dimension one) effective Lagrangian takes the canonical form 

C% ^ iTy{B[vD,B]) + F1:i{BS^,[u^,B])+DTy{BS^{u^,B}) , (2.3) 

where is the four- velocity of the baryon subject to the constraint v"^ ~ 1, Tr denotes the trace in flavor space 
and D and F are the leading axial- vector couplings, £> ~ 3/4 and F ~ 1/2. Furthermore, is the spin- vector 
and — iu^V fjUu\ where U = collects the Goldstone bosons (for more details, see U). For the calculation 
of the self-energy (mass), it sufRces to use the partial derivative 9^ instead of the chiral covariant derivative 
Df^. The dimension two chiral Lagrangian can be decomposed as (for details see QSl and JT]) 

19 

,(2) _ ^(2,br) _^ , _^(2) _^ ^(2,rc) 



with 



= 4B"^+E^^of)+4^-) (2.4) 



4b'"'^ ^ fez,Tr[B{x+,S}]+6i.Tr[B[x+,S]]+6oTr[BS]Tr[x+], (2.5) 
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hOf^ = 6iTr[5[M^, K, B]]] + &2Tr[B[u^, {u^', B}]] + h^Tv[B{u^,, {u'', B}}] 

i=l 

(&4 - mofei5)TOo + ^(^12 - »7io&i8)mo jTr[S[i; • u, [w • u,B]]] 
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+ (&5 + be- mo6i6)TOoTr[i?[t; ■u,{v ■ u, B}]] 

+ ^(67 - mobij)ma + ^(^12 - mobis)^moTT[B{v ■ u,{v ■ u,B}}] 

+ bsTT[BB]TT[u^'u^,] + bii2ie^"'°''^TT[Bu^,]vaSf}TT[u^B] 

{bg - TOofei9)rno - ^(^12 - mo6i8)mo^ Tr[SJ3]Tr[w ■ u v ■ u] 

+ 6i32ie'^''"'^Tr[Sz;„S'^[[u^,ii,],B]]+6i42ze^''"''Tr[Bi;„^^{[u^,u,],S}] , (2.6) 

with Too the octet baryon mass in the chiral hmit. Exphcit symmetry breaking embodied in the external 
source x+ ^ (with A4 the diagonal quark mass matrix) only starts at this order, collected in It is 

parameterized in terms of the LECs bo, b^ and bp. Throughout, we work in the isospin limit to,; = = to and 
thus consider four different baryon states, the nucleon doublet (N), the lambda (A), the sigma triplet (S) and the 
cascade doublet (S) (Isospin breaking corrections are discussed in ^7] and [H]). The operators with the LECs 
bi (i = 1, . . . , 19) only appear as insertions in fourth order tadpole graphs (see fTy for a detailed discussion). 
Note that in |15| the contributions from various combinations of dimension two LECs were effectively subsumed 
in one corresponding coupling since operators ^ and ~ (v ■ k)^ lead to the same contribution to the baryon 
masses. We can only do that at later stage in the calculation so as to be able to consistently work out the 
renormalization of these dimension two operators. As will become clear later, while in DR the bi are finite 
numbers, this is not the case if one employs cut-off regularization. Finally, we remark that the recoil terms 
collected in C^^g'^'^ are given in '15''. Similarly, there are further recoil corrections ~ I/'^b collected in C^j^g'^'^ ■ 
These involve no unknown parameters, their explicit form is also given in ,15, . To end this section, we give the 
fourth order terms relevant for our calculations, 

4b = WB[x+, [X+,B]]) + d2TT(B[x+, {x+, B}]} + dsTT{B{x+, {x+, B}}) + d^Tr{Bx+)Tr{x+B) 

+ d5TT{B[x+,B])Trix+)+drTTiBB)TT{x+)Tr{x+)+dsTTiBB)TTixl) . (2.7) 

We remark that we have employed the notation of jl5| to facilitate the comparison with that work and also use 
some of the LECs determined there. 

2.2 Regularization schemes 

We briefly recall the salient features of the various regularization schemes employed in calculating the baryon 
masses. Heavy baryon CHPT together with DR was used e.g. in the early papers ^1 and a fourth 

order calculation using infrared regularization (which is also based on DR to deal with the UV divergences in 
the loop graphs) was reported in JTj (for an earlier incomplete calculation, see and a recent calculation in 
the extended on-mass shell renormalization scheme was reported in JHI)- To be definite, consider the leading 
one-loop pion graph for the nucleon mass (the sunset diagram with insertions from the leading order Lagrangian 
which is of third order). In the heavy baryon approach, it is given by 

/^ = £j(0)M^ ^(0)-^/(g, (^2_,2_L)(,.fe_,,) ^ (2.8) 

with c — {D + F)'^/Fq and Fq the pseudoscalar decay constant in the chiral limit. In DR, the loop function 
J(0) is finite and can be expressed as 

J(0) = MtH^n)-^^' r r = , (2.9) 

with d the number of space-time dimensions and we have set = 4 on the right-hand-side of Eq. (|2.9|l . This 
gives the time-honored leading non-analytic contribution 

In = -i^M^.- (2.10) 

Note that in DR no power-law divergences appear and therefore loop graphs can not renormalize the baryon 
mass in the chiral limit and the dimension two LECs which leads to self-energy term ~ AI^. If we instead use 
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a three- momentum cut-off as suggested in [201, the same diagram leads to the expression 

A ti 
d\k\^ 

(P M2 -h J/2)3/2 

Note that besides the contribution ~ that is free of the cut-off, we have additional divergent and finite 
terms. The cubic divergence independent of the pion mass leads to a renormalization of the baryon mass in 
the chiral limit whereas the term ~ M^A renormalizes the dimension two LECs (the precise relations between 
the bare and the renormalized parameters will be given in the following section). Also, the finite arctan term is 
formally of higher order since it only starts to contribute at order M^. Chiral symmetry has been manifestly 
maintained by this procedure since no structures besides the ones already appearing in the effective Lagrangian 
are needed to absorb all divergences (see also [201 and the more systematic work reported in |22]). Also, the 
DR result can be formally obtained when letting the cut-off tend to infinity. This will also be discussed in more 
detail below. 




3 Formalism II: Baryon masses 

This section contains the basic formalism to calculate the baryon masses to fourth order in the chiral expansion. 
We briefly discuss the third order result and the introduction of an improvement term as proposed in 20 . 
We then proceed to present the central new results, namely the baryon masses to fourth order utilizing cut-off 
regularization. The calculation of the baryon self-energy and the corresponding mass shift at a given order in 
the chiral expansion is briefly outlined in App. ^ 



3.1 Baryon masses at third order 

The calculation of the baryon masses to third order in cut-off regularization is straightforward. Utilizing the 
loop integrals collected in App. El one obtains 



a^M^ arctan + a^M^ arctan + ct^^I^ arctan 



+ Oiq^), (3.1) 



where the state-dependent coefficients 7^'^ and can be found e.g. in As announced, the baryon 
mass and the the dimension two couplings are renormalized as symbolized by the superscript (r). The precise 
renormalization takes the form 

(r3) f \ ^3 .(r3) , / 5DF \ ^ 



(.,3) , / 13I?2 + 9F2x ( -D^ + -iF^ \ 

bo = bo~(^^^^)A, =bn-[ 3^^,^, jA. (3.2) 



It is instructive to expand the arctan(Mp/A) {P = {tt, 77}) contributions 

= mi'-^+^Eb^^^+j^bP-2b^;\M'^ + 2M],)-^[alM'^ + a^A^ 



127r2i^o " 1 A 3 V A 



lEb%^ + iSbP - 2bt\Ml + 2Ml) - [alMl + a^M^ + a^M^] + 0(g4) ,(3.3) 
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where the last hne corresponds to the DR resuh. As stated earher, the additional contributions scale with inverse 
powers of the cut-ofF and thus vanish when A —^ oo (which formally corresponds to DR). As already stressed in 
j20| (and noted by others), this third order representation is not sufficiently accurate to make connection with 
lattice results if one is not very close to the physical value of the GB masses. Therefore, one should perform 
a fourth order calculation or at least add an improvement term that is formally of fourth order but should be 
elevated to third order. We first discuss brieffy this latter possibility before turning to the fuU-ffedged fourth 
order calculation. 



3.2 Improvement term 

As noticed in (20|, the third order result for the baryon masses shows a very strong cut-off dependence when 
the pion mass is increased above its physical value, see the left panel in Fig. ^ Only when one has a plateau 
below the chiral symmetry breaking scale = AttFt^ ~ 1.2 GeV, one has the required cut-off independence. 
For pion masses above 300 MeV, this plateau vanishes. To obtain a better stability against cut-off variations, it 
was therefore proposed in (20] to promote the fourth order operator eiM^NN to the third order and to cancel 
the leading cut-off dependence in Eq. H3.3|l by a proper adjustment of the LEG ei, ei = e^" — coeff/A, where 
the coefficient can be read off from the leading term of the expansion of the arctan function. For the SU(3) 
calculation performed here, the situation is a bit more complicated. In fact, the corresponding improvement 
term for the baryon B consists of three contributions (we use the notation of 

el%Mt + el$MlMl + e^^f Mj, , (3.4) 

where the coefficients e^'^ (P, Q = {vr, K}) are linear combinations of the fourth order LEGs di defined in 
Eq. H2.7|l (the precise relation can be found again in ^H])- Throughout, we use the Gell-Mann-Okubo relation 
to express the M^-term by the pion and the kaon masses, SM^ = 4M|- - Ml. To ehminate the leading 1/A 
dependences, the LEGs di have to take the form 



iJ^-3F^ DF _ D^~iF^ _ -D^ + iF' 

imF_ , .fin , 35^!+27F! -rffin , 17^!+9F! .3 .^ 



3.3 Baryon masses at fourth order 

To fourth order in the chiral expansion, the octet baryon masses can be written as (when employing GR) 



(2) 
B 



,(3) 



,(4) 



4^ 



Ami'^ + /b,3(A) + Amf + /s,4(A) 



(3.6) 



where in the second line we have split the mass shift Sra'j^ (i = 3,4) into cut-off independent and an explicitly 
cut-off dependent piece. This is done to facilitate the comparison with the results obtained in DR. The various 
pieces take the form 



5m 



(2) 



A? 



,(3) 



/b,3(A) = 



247ri^2 
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a^A/^ arctan + a^M% arctan — h al,M^ arctan 



A 



(4) 



Am 



€^^§MM + 4§,tMImI In 



-4;lj.MM\nil + RT) 




A 



A 



l.B.T 



AM 1 - i?T 



1 {Mt^ 

T 



+P^,b.tM''pA:' (1 - i?T) + P^'3,tM''pM-'q 



(3.7) 
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with Rt = (1 + M|,/A^)^/^. Here, we have P,Q,T — {tt, K, 77} and the summation convention for these indices 
is understood. To arrive at these results, we made use of the loop integrals collected in Add. IbI The resulting 
A-dependent terms are separated into a contribution that contains all the power and logarithmic divergences 
and another one that is finite in the limit that A —^ 00. Only these latter terms are displayed here. For the 
nucleon, the e-coefficients and the /3-coefHcients are collected in App. Note also that the fourth order LECs 
di have the form displayed in Eq. (|3.5(l so that the leading cut-off dependence is canceled. In the second order 
term, we have used the leading terms in the quark mass expansion of the pion and the kaon mass, denoted 
Mo,p, because the difference to the physical masses appears in the fourth order of the baryon mass expansions. 
For consistency, we have recalculated the Goldstone boson masses in cut-off regularization (this was not done 
e.g. in the explicit formulae are given in Add. IO All appearing polynomial and logarithmic divergences 

are absorbed in a consistent redefinition of the bare parameters. The renormalization of the chiral limit mass 
and of the LECs takes the form 

(rA) _ (r,3) , mp / 612 _ 864 _ Tby Sbismp Ibnmo bismp \ 4 



rrin = rrir, H t, — ^ Oq H h Oi qmn A 

° ° t:^F^ V 8 4 12 ^ 4 12 8 V 



^(r,4) _ f,(r-,3) 

-bb^mp — bb^mp -\- Sfoigmg ) A^, 



b'p' ' = b'p'> - ^ ( - IO62 + 145f + mpD' + 206i)DF + 186^^^' 



M) _ ^{r.z) / f - i8bi - 26b3 - 48^8 + IS/^g + 486o - I'obuD^ - mbpDF - ISbcf ^ 



K - b 

Oq — Oq 



9&4TO0 ~ IS&vTOo ~ 246gmo + 96i5TOq 136i7mQ -1- 246igmQ j A , 

r,4) _ ;,(r-,3) 

ISbijF^ - 36i2mo - 964TO0 - bjmo + 96i5mg -I- 6i7mg + Sfoigmg ) A^, 



bo"' = ^D^' - \ F2 ( - - 263 + 145c + 266ci?2 ^ ^qj^^dF 



4"'*^ = rfi - \ {- D'' + iF^ + l2bimp-Ab:imo-UbDmp 

11527r^FQ"^mg V 

—69bDD^mQ — l(S2bpDFmQ — SlbjjF^mo + SbiUiQ — feymg — 36i5TOq + 617^0^ In , 

/ Too 

4"'*^ = d2-—4^^ ( - 6Z?F - 1262mg - 46Fmo - eObj^i^^mo 

YOSTT^i'g TOg V 

A 



-UObDDFmo - lOSbpF^niQ - ib^ml - 3beml + 3bieml ) In ■ 

/ Too 

4''"''^ = d3 - {~iD^ + 9F^ + 365iTOo + 463TO0 - 246dtoo - ISboD^mo 

YOSTT^i'gTOg 

— lOSbpDFrriQ — bAbjjF^mQ + 36i2TOq + 9fo4mQ + fermg — 96i5TOg — 6i7TOg — S&isTOg) In , 

y Too 



4"'*^ = di - ^^^^ ^ „ (9D^ ~ 27F^ - IO861TO0 + 453TO0 + 445dTOo 

11527r^FQmg V 

-h288&Di:'^TOo - 6&12TO0 - 2764TO0 -I- 67TO0 -I- 27bi5ml - bnml + Qbiaml] In — , 

/ TOg 

4''''*^ = 4 - 1^2 (2613^ + 5262TO0 - 446ftoo + 1365^1^ + meml - 136i6TOg) In , 

dp^^ = d7- ^^^^ ^ „ (35D^ + 27F^ + IO8&1TO0 + UObamo + 264bsmo 

69127r^Fg"'mg v 

-132&£,TOo - 2646gTOg + lAAboD'^TTiQ + 27^4 to^ + 35bTml + 66bgml 

— 276i5TOg — 356i7TOg — GGbigTOg") In , 

/ Too 

4'''''^ = ds- l^n (l7D^ + 9F'^ + 366iTOo + 6863^0 + 12058mo 



23047r2Fn2 



qTOo 
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~28bDmo - I2O60TO0 + leSboD'^mn + 4:32bFDFmo + 216bDF'^mo 

+9b4ml + nbTinl + SOfogmg - 96i5mo - nbnml - SO&igmo) In — . (3.8) 

/ Too 

Clearly, higher order powers in the cut-off A appear as compared to the third order calculation. Note also the 
appearance of a new scale in the logarithmically divergent terms. As can be seen from App. ^ the integrals 
h — /g and ai — contain terms proportional to ln(Mp/A). To properly absorbs these divergences, one uses 
ln(Afp/A) = \tv{Mp/v) — ln(A/i/), with v the new scale. From here on, we set v = TOq. As a check, it can be 
shown that for A — > cx) our fourth order results agree with the ones obtained in DR in |15| if one sets fj, = toq, 
with /i the scale of DR. For this comparison to work, one also has to account for the fact that in jTH] some terms 
~ MpMg were absorbed in a redefinition of the di. This concludes the necessary formalism, we now turn to 
the numerical analysis. 

4 Results and discussion 

Before presenting results, we must fix parameters. In the meson sector, we use standard values of the LECs Li 
at the scale Mp-. = -0.3, L5 = 1.4, Lg = -0.2, L7 = -0.4 and Lg = 0.9 (all in units of 10"^). These are run 
to the scale A = TOq using the standard one-loop /3-functions. Throughout, we set Fq = 100 MeV, which is an 
average value of the physical values of F^, Fx and F^. For the leading baryon axial couplings we use I? = 0.75 
and F — 0.5. We have checked that varying these parameters within phenomenological bounds does not alter 
our conclusions. 

4.1 Fixing the low-energy constants 

For the dimension two LECs from the meson-baryon Lagrangian we use the central values of |15) . these are 
collected in Tab.^ We have not varied these LECs since such modifications can effectively be done by changing 
the fourth order LECs di within reasonable bounds. Next we consider the determination of the fourth order 



60 


ho 


hp 


hi 


h2 


63 




-0.606 


0.079 


-0.316 


-0.004 


-0.187 


0.018 


-0.109 



Table 1: Values of the LECs h in GeV"^ taken from [El. 



LECs di. As noted before, the improvement term for each baryon consists of three pieces. To determine these, 
we have only considered the nucleon, for two reason. First, there are more lattice results for this particle than 
for the others and second, it also facilitates the direct comparison with the SU(2) results of [201 • For the nucleon, 
the coefficients appearing in Eq. (|3.4|) are related to the LECs di via 

^Vn = -4(4^1 + 2^5 + ^7 + 34) , ej^ = 8(4^1-2^2-4-2^7 + 2^8) , eff = -16(^-4 + 4-4 + 4 + 4) • 

(4.1) 

We have now varied the values of the di under the following constraints: We require that mM{MT^) passes through 
the physical value ruN — 940 MeV for the physical pion mass ilf^r = 140 MeV and similarly for mN{MK) at 
Mk — 494 MeV. Furthermore, we only allow for variations of 5di — ±0.1GeV~^ from the central values of 
[T^ (this is in fact the largest magnitude of any of these LECs). Under these restrictions, we have tried to 
describe the trend of the earlier MILC data with the third order improved formula (note that the more recent 
data 01 were not used in the fit for reasons discussed below). Note that we have reconstructed these data from 
table IX (VII) of Ref. 0(131) using the scale parameter ri — 0.35(0.317) fm. The resulting set of values for the 
LECs 4 is denoted as the "optimal set" from here on. In Tab. |21we have collected the values of the 4 from 
and for the optimal set. In fact, the values of the di given in that table refer to the basis used in TS' as 
indicated by the superscript "BM" . These differ by some small finite shifts from the one used in CR (we refrain 
from giving the explicit formulae here) . Note that one can not exactly reproduce these lattice data as shown by 
the solid line in Fig.|21 but the nucleon mass now increases with growing pion mass as demanded by all existing 
lattice results. Also, at fourth order there are other contributions which arc not captured by the improvement 
term, see the discussion below. We also note that we have not restricted the 4 such that the CMC relation 
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3mA + ms = 2(r7ijv + ?tih) is fulfilled so as to have a better handle on the theoretical uncertainty. Including the 
improvement term with these values of the LECs leads indeed to a very reduced cut-off dependence as shown 
in the right panel of Fig. ^ Note that the treatment of the improvement term is more tricky in SU(3) than in 
the SU(2) calculation since one has to balance three different terms as opposed to fixing one in the two-flavor 
case. 







dT 




dr 


jBM 
"5 




jBM 
"8 


BM 15 


0.008 


0.035 


0.069 


-0.077 


-0.05 


-0.018 


-0.103 


opt. 


-0.043 


-0.066 


-0.031 


-0.077 


-0.15 


-0.118 


-0.2 



Table 2: Values of the LECs di in GcV in the basis used in as indicated by the superscript "BM". "Opt. 
denotes the optimal set as described in the text. 



4.2 Nucleon mass and pion-nucleon sigma term 

We now consider the nucleon mass as function of the pion and the kaon mass in DR and CR. We have calculated 
TTiAT at third, improved third and fourth order, see Fig.|2for DR. In what follows, we mostly focus on the results 
obtained at fourth order. The trends when going from third to improved third to fourth order for the pion 
mass dependence of the nucleon mass are very similar to the SU(2) case discussed in big detail in [201 ■ We see 
that with the optimal set of the di as given in Tab. El one obtains a rather accurate description of the earlier 
and most of the more recent MILC data [H] for pion masses below 600 MeV, cf. the dot-dashed line in Fig. |21 
Note, however, that the two lowest pion mass points of the recent MILC data |3j do not quite fit into the trend 
of our extrapolation function if one insists that for the physical pion mass the curve runs through the physical 
value oi niM- More low mass pion data and/or a more sophisticated treatment of finite size/ volume effects 
are needed to resolve this problem. If one were to use the di determined in [151, one already deviates sizeably 
from the trend of the MILC data for pion masses starting at about 500 MeV (dotted line in Fig. [SJ . The same 
can be seen for the fourth order calculation based on CR utifizing A — 1 GeV in Fig.|21(left panel). To get a 
better idea about the uncertainty when going to higher pion masses, we have also performed calculations with 
three other sets, namely setting all di ~ 0.2/0/ — 0.2 GeV""^, corresponding to the long-/medium-/short-dashed 
lines in that figure. This clearly overestimates the theoretical uncertainty since some of the di are correlated 
parameters. Still it is safe to say that for pion masses below 550 MeV the theoretical error is moderately small. 
These results for the pion mass dependence of mjv as well as for its cut-off dependence at a given pion mass 
are very similar to the results of two- flavor study reported in [201 • In the right panel of Fig. O we show the 
kaon mass dependence for the same variety of choices for the di. Since we enforce that rriN takes its physical 
value for Mk — 494 MeV, the resulting kaon mass dependence is much flatter than the pion mass dependence 
with decreasing meson masses. In the left (right) panel of Fig. 01 we show the cutoff dependence of for 
various values of the pion (kaon) mass. For pion masses up to 450 MeV, one has a nice plateau below the 
scale of chiral symmetry breaking but not any more for M^r = 600 MeV. For the kaon mass dependence, the 
situation is somewhat different due to the much larger meson mass. Here, we still have a reasonable plateau at 
Mk — 600 MeV. These observations are consistent with our earlier observations that chiral extrapolations in 
based on the fourth order CHPT representation can be applied for masses up to 550 MeV, a result which 
is consistent with the one found for the SU(2) calculation in [201 ■ -F'^'" kaon mass dependence, one can even 
go to somewhat higher meson masses. 

It is also interesting to study the range of values found for the octet chiral limit mass mp at the various orders 
and employing the different regularization schemes and values of the LECs di . One observes that mo increases 
with increasing cut-off, that means in DR its value is above the one in CR when one chooses A = 1 GeV. 
Insisting that mjv takes its physical value at the physical value of M^r and Mk when studying the pion and the 
kaon mass dependence, respectively, we find 

710 MeV < mo < 1070 MeV , (4.2) 

which is consistent with expectations and also with the findings in |15j (note that in that paper a different 
method was used to estimate the theoretical uncertainty, which we consider less reliable than the one used 
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here). The range given m Eq. 14. 2|) of course includes the SU(2) value of about 880 MeV [201 ■ Also, we note 
again that the MILC data are obtained using staggered fermions, so strictly speaking one should use "staggered 
fermion CHPT" . However, we believe that this will not significantly alter the trends discussed here. 

Another quantity of interest is the pion-nucleon sigma term, 

a^^(0)^ni{N\uu + dd\N)^m^^M^^ . (4.3) 

It can be extracted directly from the the slope of mN{M^) at M^r — 0. For the optimal set and the LECs from 
|15| . we obtain the range (considering DR and CR with A = 1 GeV) 

cr^^ (0) = 50.7... 53. 7 MeV . (4.4) 

These numbers are consistent with the results of 15 (as they should) and the study of SU(2) lattice data in 
|23|. (T^AT = 49±3MeV. The resulting strangeness fraction can be obtained from (T7rjv(0) = ao/{l — y) and using 
(To = 37 MeV from ^| (which is consistent with the pioneering work in [21], cto = 35 ± 5 MeV). This leads to 
y = 0.27 . . . 0.31, which is again consistent with but somewhat on the large side. 



4.3 Hyperon masses 

We now consider the octet members with strangeness. As noted before, when fixing the coefficients in the 
improvement term, we have not insisted to recover the Gell-Mann-Okubo relation, thus some of the masses are 
somewhat off their empirical values. In Tab.|3|we collect the resulting values for the improved third and fourth 
order. While the E mass is well reproduced, the A and S masses come out by about 10 — 15% too high. To 
get a handle on the theoretical accuracy, we also use the values for the di from ^H], in that case all masses are 
exactly reproduced. 



order / 


imp. 3rd 


fourth 


fourth 


exp. 


baryon 


CR 


CR 


DR 




A 


1115 


1304 


1243 


1116 




1101 


1194 


1167 


1193 




1222 


1532 


1437 


1315 



Table 3: Baryon masses in MeV in DR and CR with A = 1 GeV for different orders. For the experimental 
numbers, we haven taken the masses of the neutral particles. 



The corresponding pion and kaon mass dependences for the A, the E and the S are shown in Figs. |S1 [HI and 
[TJ respectively. The solid/dashed lines refer to the optimal set of the LECs/to the LECs from ^Jl- We note 
in particular that the pion mass dependence for the S is much flatter as one would expect from the MILC 
data. This is not unexpected - the S only contains one valence light quark and should thus be less sensitive to 
variations in the pion mass. Clearly, one could improve this description by fitting directly to these particles. 
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A Baryon masses 

Here, we collect some formalism to calculate the baryon masses from the baryon self-energies. Consider the 
heavy baryon approach. The baryon four-momentum is — niQV^ + r^, with t;^ the four- velocity subject to 
the constraint = 1 and is the (small) residual four- momentum, uj = v ■ r <^ toq. The baryon self-energy 
E(aJ, r) has the chiral expansion 

^{uj,r) ^Y.^^\uj,r) ^Y.'^''\uj,r) + i:'^^\uj,r) ^ . . . , (A.l) 
and the corresponding baryon mass shift is then given by (with Sitib = Sirig^ + Sra"^^ + 5rri"g +...): 

<5to^^' I](2)(^^0,r) + -^, (A.2) 



(2) 
B 








(3) 
B 


= I](3)(c^^0,r) + ^S(2)(^ 
auj 


= 0,r) 


(4) 
B 


2rno 


+ A: 


doj 




+ ^S(2)(^^0,r)5mg) + 


1 92 




29^ 



,(2) _ 

2mo 



6m^g> = I](3)(cj^0,r) + ^S(2)(^^0,r)((57r4'^- — ), (A.3) 



(A.4) 



(2)^2 c (2) r 



tuq Arrvk 



Prom this, one obtains the pertinent representations of the baryon masses. In the following, we only discuss 
the nucleon mass. More precisely, we now give the corresponding non-vanishing prefactors for the nucleon (we 
refrain from giving the coefficients of the other octet members). At third order, cf. Eqs. one has the 

standard values 



= -AD + F)\ a% = -{5D'-6DF + 9F') , al^-{D-3Ff. (A.5) 



At fourth order, see Eq. H3.7|l . wc have a cut-off independent and a cut-off dependent contribution to the nucleon 
mass shift. The coefficients of the A-independent term read 

el^^ = -16di - 8^5 - 4d7 - 12^8, 

= 32di - 16^2 - 84 - 16d7 -I- ledg, 
e^^ = -I6di + 16^2 - 16^3 + 16(^5 - 16^7 - 164, 

1 / 3Z)^ (}DF 

. = TTTT-Er^ - - 12^2 - I263 - 2468 + 125d + 126f + 246o 

' ' 2(47ri^o) V ^0 mo 

3^2 

3647710 — 3657^0 — 3667710 — 3677^0 — 6697770 + 3615TOQ 

7770 

.2 I 01, ™2 



-|-36i6777q -t- 36177^0 -I- 66l9777g 

4nk = , ^ ,J -52bDD^ + 60bFD^ + 120bDDF -72bFDF -36bDF^ + 108bFF^ 
' ' 3(47ri^o)'^ V 

<^2NK = 777— ^77^1 - 366i-f 1262 -3663 - 4868 -H366d-126f + 486o + 526i5i:'2 



3(47r^^o)2 



5r)2 

-eObpD^ - UOboDF + l^bpDF + SGboF^ - lOSbpF^ 



7770 

6DF 

36i277io — 9647^0 + 36570,0 -I- 3667^0 — 9677770 

7770 mo 

-12697^0 -I- 96i5rr7Q — 36i67r7g + 9bnmQ + Sbism^ + 12bigmQ j . 
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""Tnt, = ^TT-^-^TTo f - 366i + 1262 - 463 - 2468 + 366j3 - 606f + 246o - — 
' ' ' 54(47rFo)^ \ mo 



6DF 



mo mo 

964mo + S^smo + S^emo — brmo — Sftgmo + dbi^m?, 

mo 



-36i6TOo + bnml + Qbigml 
1 / „. „„. „. „. 8D' 



,2 



28861 - 9662 + 3263 + 19268 - 2406d + 3366f - 1926o + 



16D 



2 



54(47rJ^o)n . ---o , ^^-^o , , 

48£)F 72F2 

1 h 7264mo - 2465mo - 2466mo + 867mo 

mo mo 

+4869mo - 726i5mo + 246i6mo - 86i7mo - 486i9mo 

^2N n = ( ' 5766i + 19262 - 6463 - 38468 + 3846^ - 3846ir + 3846o - 

' ' ' 54(47rFo)2 V "^0 

96DF 144F^ 

H 14464mo + 4865mo + 4866mo - 1667mo 

mo mo 

-9669mo + 1446i5mo - 486i6mo + 166i7mo + 966i9mo^ . (A.6) 

The corresponding coefficients of tfie A-dependent term are 

A.jv.TT = J^^^ - 364 - 365 - 366 - 367 - 669 + mo(36i5 + 36i6 + 36i7 + 6619)^ , 

Pi,N,K = ( - 612 - 364 + 65 + 66 - 367 - 469 + mo(36i5 - 616 + 36i7 + 618 + 4619) 

8(7rFo)^ V 



Pi,N,n = ^(^^p^y ( - 964 + 365 + 366 - 67 - 669 + mo(96i5 - 36i6 + 617 + 6619) ) , 



/52,jv,7r = FN2 ( - 661 - 662 - 663 - 1268 + 66d + Qbp + 126o - 364mo - 365mo 



-Sbemo - Sbrmo - Gbgirio + 36i5mo + 36i6mo + 36i7mo + 66i9mo ) , 



1^2 NK = ^,0 -52bDD'^ + 60bFD'^ + 120bDDF -72bFDF -SeboF^ + 108bFF^ 

' ' 9(47ri^o) V 

P2N K = V ^■2 ( - IO861 + 3662 - IO863 - 14468 + 1086d - 366f + 1446o + 526d-D^ 
' ' 9(47rFo)^ V 

-eObpD^ - nOboDF + 72bFDF + SGboF^ - lOSbpF^ - 186i2mo 
— 5464mo + 1865mo + 1866mo — 5467mo — 72bgmo + 546i5mo 

-186i6mg + Mbnml + 186i8mg + 726i9mg^ , 

/32,jv,T, = Q(^4^p^Y (l^^i ~ + 263 + 1268 - 186r, + 306f - 126o + 964mo - 365mo 
— 366mo + 67mo + 669mo — 96i5mQ + 36i6mo — bnniQ — 66i9mQ^ , 

1^2 N n = V x2 I - 726i + 2462 - 863 - 4868 + 486d - 486f + 486o - 3664mo 
' ' ' 9(47rFo)^ V 

+1265mo + 126f;mo — 467mo — 24:bgmo + 366i5mQ — 126i6mQ 
+46i7mg + 246i9m^ ), 
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1 / 3I?2 

PTn . = - - 12&2 - 1263 - 2468 + 126d + 126f + 246o 

' ' 2(47rFo)^ V 7710 

364777,0 — 3657710 — 3657770 — 3677770 

7770 "lo 

—6697770 + 36i5777o + 36i6777o + 36i7777o + 66197770^ , 

1 

3(4^o)2 
1 

3(4^o)' 



f^l%,K = ^^ ( - 526i5L'2 + mhp + 1206d DF - l^hpBF - 366d-F^ + 1086^^^^ 

Pt^^K = TTTj^TTT^ ( - 366i + 1262 - 3663 - 4868 + 366d - 126f + 486o + 52bDD'^ 



^eObpD^ ~ UOboDF + nhpDF + 366Di^^ - impF"^ 



5D^ 



Q7T-K 

P: 



mo 

6DF 9F'^ 

H 36127710 — 9647710 + 3657710 + 3657770 

7770 "^0 

— 9677770 — 12697770 + 96i 577tQ — 36ig777Q + 96i7771q 
+ 36l8777,o^ + 126l9777o^ , 

263 468 26£, 106f 46o 



(47rFo)2 V39 27 9 3 9 9 54777o 

647770 br^rriQ barrio bjmo 



-+ 


9 : 


DF 




97770 


67710 


69777-0 


6i5777j 


9 


6 



6 18 18 54 

2 h._^2 u^^2 



18 54 9 

_ 1 /^166i I662 , I663 ^ 3268 4 06d , 566f 326o , 4D^ 



(47rFo)2 V3 9 27 9 9 9 9 27777o 

8DF 4:F^ 4647770 4657710 4657770 4677710 
~ 97710 ^ 3^ ^ ~3 9 9 ^ 27 

8697770 ^ 46l5777§ 46ig777^ ^ 46l 7777g ^ 86l 9777|\ 

9 3 9 27 9/' 

_ 1 / 326i 3262 3263 6468 6 46d 646^ 646o SD"^ 
^3,N,^ - (4^_p^)2 l^"~ + ~9 27 9~^ ~9 9~ ^ ~9~ " 27^ 

16DF 8F2 8647770 8657770 8657770 8677770 
97710 ~ 3^ 3~ ^ 9 ^ ~~9 27~ 
I6697770 8615777^ 8615777^ 86i7?77g 16619777^ ^ 

Finally, we note that we use the values of the bi as given in Table Q] and set all other dimension two LECs to 
zero. 



B Loop integrals in cut-ofF regularization 

In this appendix, we collect all integrals needed in the calculation of the baryon self-energy to fourth order 
utilizing cut-off regularization. Here, M is a generic symbol for the propagating Goldstone boson and we give 
the explicit representations only for the relevant case M > \u}\: 

h(M,Lo) = ( S^-. -o ^ -{S-kf 




— - (M^ - cj^) A + ( Af ^ - cj^ 1 TT - ( - uj^ 1 arctan ^ 
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wA^ r M2 / 3 \ /A 



+ \ 1 + 



A2 



h{M,0) = - 



-(M^ - arctan ( ^/ - 1^1 + ^ 



--M^A + M3--M3arctan - 



(B.l) 



h{M,0) 



h{M,oj) 
h{M,0) 

h{M,0) 
h{M,oj) 
h{M,Q) 



:{S-kf{v.k) 



h{M,0) = 



—I 



(27r)4 - + iO+ uj-v ■k + iO+ 



A 
M 



A2 
M2 



IStt^I- 



A*Wl + :f---M^A^/l + 



A2 8 



A2 



it d 
(27r)4 - M2 + zO+ (a; - • fc + tO+)2 ' = "'dJ^^^^''^) 



^{2 

167r2 I 2 



M2 

A2" 



A^/l + ^+3a;A- 



A2 



M2 



-3.Vm^. + (3.2 _ ^M2) m + f7§,) 

+3w\/ — a;2 arctan f 



+3a;-\/ - arctan ( 



V A 

M2~ 



1 fA2 f M2" 



\ V 

1 



A2 / ' 



1 + 



2 



y (27r 



TT-n5 ^4; — ^7 ^ ^rrT^(5' • kf{v ■ k) = -i-^l2(M,uj) 



i7i(M,w) +a;73(M,a;) 

i/i(M,0) , 



(27r)4 /c2 _ ^2 _^ ^0+ (w - w • fc + zO+)2 
iujIi{M,uj) -il2{M,Lv) +uj^h{M,u)) , 
i/2(M,0) , 

d^fc Z z2 

(27r)4 A:2 _ 71^2 + ^0+ (w - t; • fc + iO+)2 



(5-fc)2(t;. A;)2 



(S -kyk^ = M^h{M,Lu) , 



= M2/3(M,0) . 



+ LOh (M, w) 



(B.2) 
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Similarly, for the calculation of the baryon tadpole and the meson masses in CR we need the following integrals: 

1 f . I M2 

In 



a2{M) 

Q!3(M) 



2(27r)2 

J {27T 



A2 



mo J 



-In 



mo J 



i 



-k^ - 



M^ai{M) 
1 1 



rAS 1 



(27r)4 F - M2 + iO+ ° (27r)2 4 

C Meson masses in cut-off regularization 



Af2 

A2" 



-M'^ai{M) . 



Here, we collect the formualae for the meson masses to fourth order in CR. The pertinent diagrams are tree 
graphs with one insertion from C^? and C^f^ and tadpoles with exactly one insertion from . We have 



Ml = M^^ + 5Mi^^ 

TT U.TT ' TT 



M4 ^^M^ 



+ 



+ 



WF^w'^ mo 48Fgn^ mo 
M.4 ^ 



+ 



Fi 



Ml 



lnll + Wl+(^ 



F^ 

" "-In 1 



48Fo^7r2 



1 



2 I 16J;f7r2 



A 



+ 



16F27r2 



M2 



48Fo^7r2 



AM1-W1 + 



A 



M|M2 
48Fo27r2 



(C.l) 



MtM^ Mr, MtMl , lr\ (r\^ Mi , (r\ (r-\ (r-\ ^r1^ 








2 



_ M2.M2 

24F2^2-1^-Y'"V a; i^^/7-7KT^l24FoV^^ V' y^\AjJ 2AF^-,^ 



(C.2) 



= M„%-^ln^ + J^ln^+f-^^ + H^-^)ln^ 
16^1)2^2 ^ 6^2jr2 TOO V 432F27r2 108Fo27r2 27F^n^J mo 

Mi, 64^1''^ 128L1'^^ 128l1''^ 128L^.''^ 128l1''\ 

J £lf ± 2 1 2 1 L 1 2_^ 

^ F^^ 3 9 ^ 3 ^ 3 ^ 3 ^ 

M|M| 16L^ 64L^''> 32L^;'> 256L^;^ ^^SL^sK 
^Fo 3^9^3 3 3^ 

SL^) S^r^ 164''^ , 128ir , 1 fi, (.)^ 
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WF^TT^ \ y \ A J GF^TT^ "' \ ^ ' V ^ ' V A 



7M4 llMlMl AMI, ^^^(^ , , (M,y 



432i^2^2 i08i^2^2 27F^TT^' \ ]] \ A 



M„\^] WF^TT^ \ \ \ A J J UF^TT 



2^2 



(C.3) 



1 Afi , ^M_^^^2f^_ L^fMlV\ 4M| , llAf^Mi 7M, 



M„^2P 9i^27r2 144i^2^2^ \^ y \A J J 27F^n'^ lO&F^n'^ 4:S2F^tt^ 



As required, in the limit A ^ oo we recover the standard DR resuh 25 . The polynomial and logarithmic 
divergences in the cut-ofF are taken care of by the following renormalization (note again that e.g. Bq is not 
renormalized in DR): 

° " 247r2i^2 " ' ^ ^ 3 3 23047r2 mo 



4''^ - 2L^p = L^-2L^ + -^\n — , - 2L^^^ ^ - 2L^ ~ — ^ In A . (c.4) 

Odtt^ mo Syovt-^ mo 



Here, Bq connects the leading terms in the chiral expansion of the Goldstone boson masses with the quark 

masses. 

Ml, = 2B^;'^m , MIk = B^^^ (m + m,) , MJ^„ = ^B^''^ (m + 2m3) , (C.5) 
with m the average light quark mass. 
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Figures 



2 3 
A [GeV] 




Figure 1: Cut-off dependence of the nucleon mass for various pion masses with the kaon mass fixed. 
Sohd/dashed/dot-dashed/dotted line: AU = 140/300/450/600 MeV. Left panel: Third order calculation. Right 
panel: Third order calculation with the improvement term. 




[GeV] 



Figure 2: Nucleon mass in DR at third (dashed), improved third (solid) and fourth (dot-dashed) order, respec- 
tively. The dotted line represents the fourth order calculation from ^Sl- The three flavor data are from the 
MILC collaboration (boxes from ||2j and stars from ^ ). The filled circle gives the value of the physical nucleon 
mass at the physical value of Mtt. 
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Figure 3: Left panel: Pion mass dependence of the nuclcon mass for various sets of the LECs di as explained 
in the text. Right panel: Kaon mass dependence. 




2 3 
A [GeV] 




2 3 
A [GeV] 



Figure 4: Left panel: Cut-off dependence of the nucleon mass for various pion masses with the kaon mass fixed. 
Solid/dashed/dot-dashed/dotted line: — 140/300/450/600 MeV. Right panel: Kaon mass dependence for 
fixed pion mass. Sohd/dashed/dot-dashed line: Mk = 494/600/700 MeV. 
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[GeV] 




Figure 5: Left panel: Pion mass dependence of the A mass for various sets of the LECs di as explained in the 
text. Right panel: Kaon mass dependence. 




Figure 6: Left panel: Pion mass dependence of the S mass for various sets of the LECs di as explained in the 
text. Right panel: Kaon mass dependence. 





Figure 7: Left panel: Pion mass dependence of the S mass for various sets of the LECs di as explained in the 
text. Right panel: Kaon mass dependence. 
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